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ABSTRACT 

In this work, we derive the stellar initial mass function (IMF) from the superposition 
of mass distributions of dense cores, generated through gravoturbulent fragmentation 
of unstable clumps in molecular clouds (MCs) and growing through competitive ac- 
cretion. MCs are formed by the turbulent cascade in the interstellar medium at scales 
L from 100 down to ~ 0.1 pc. Their internal turbulence is essentially supersonic and 
creates clumps with a lognormal distribution of densities n. Our model is based on 
the assumption of a power-law relationship between clump mass and clump density: 
n oc m x , where a; is a scale- free parameter. Gravitationally unstable clumps are as- 
sumed to undergo isothermal fragmentation and produce protostellar cores with a 
lognormal mass distribution, centred around the clump Jeans mass. Masses of indi- 
vidual cores are then assumed to grow further through competitive accretion until 
the rest of the gas within the clump is being exhausted. The observed IMF is best 
reproduced for a choice of x = 0.25, for a characteristic star formation timescale of 
~ 5 Myr, and for a low star formation efficiency of ~ 10 %. 
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1 INTRODUCTION 

The origin of the initial mass function (IMF) is a long- 
standing issue in modern astrophysics. It has numerous im- 
plications: from cosmology (e.g. cosmic reionization and for- 
mation of first galaxies), through studies of galactic struc- 
ture and evolution, down to formation of planets and plan- 
etary systems. As such, an explanation of the IMF is one 
of the key goals of star formation (SF) research. It should 
reflect not only the variety of initial conditions in SF sites 
but implement complex and intertwined physical processes 
like gravitational collapse and fragmentation, turbulent mo- 
tions, shock waves, accretion and protostellar outflows. | 

The extensive photometric and spectroscopic studies of 
the Milky Way and some nearby galaxies over the last two 
decades led to reliable determination of the IMF in the range 
from 0.02 - 0.07 M Q (brown dwarfs, BDs) up to ~ 100 M Q . 
The high-mass part of the IMF (> 1.0 Mq) is nowadays 
firmly established to be a power-law function dN/dlgM oc 
M r . The univ ersality of its 'Salpeter slope' T ~ —1.35 
ijSalpeterl Il955l ) is confirmed by i nfrared observations of 
the A rches cluster |Kim et al.ll2006l : [Pib. Kim fc Shadmehril 
although r could vary (±0.5) in some regions of active 



★ E-mail: eircne@phys.uni-sofia.bg 



star formation (see Elmegreen 2009 for a review). The low- 
and intermediate- mass IMF (0.08 < M < 1 Mq) is a plateau 
of remarkable unifo rmity under various environmental c on- 
ditions in the ISM (|Elmegreen. Klessen fc Wilson! [20ul h It 
could be rep resented with a power-law fit with much shal- 
lower slope jKroupal 1200 il l or with a lognormal function 
(|Chabriej|2003l '). Significant uncertainties remain regarding 
the B P range of the IMF. For instance, iThies fc Kroupal 
|2007l ) demonstrate sensible discontinuity at ~ 0.08 Mq in 
several young Galactic clusters. 

Efforts dedicated both in analytical and semi-analytical 
studies and in numerical simulations led to obvious ad- 
vance in reproducing the observational IMF in some Milky 
Way clusters but the picture is far from being complete 
(cf. Bonnell, Larson & Zinnecker 2007 for review). One 
current theory describes the S F process as being con- 
trolled by supersonic turbulence dPadoan fc Nordlimd 2002; 
Mac Low fc Klessen! I200I iBallesteros-Paredes et all 120061 : 



Hennebelle fc Chabrierl |200S| , 120091 ). This process is often 
referred to as gravoturbulent fragmentation. The physics 
of turbulence is still quite poorly understood, because of 
the great mathematical complexity of the fluid equations 
(see e.g. Lesieur 1997) and the issue of how turbulence is 
driven remains under debate (|Mac Low fc Klessen! 120041 ) . 
Numerical simulations show that although supersonic tur- 
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bulence can provide global support, it produces density en- 
hancements in molecular clouds (MCs) that allow for lo- 
cal collapse - gravitationally bound clumps of a few hun- 
dred solar masses contract and fragment with forma t ion of 
compact protostellar cor es dKlessen fc Burkertl Eoool . l200ll ; 



iBate. Bonnell fe Brommll2003l ; I Clark et al .1 120051 ) 



If such forming clusters are gas-rich, competitive accre- 
tion is another crucial factor for the origin of the IMF - it 
determines the final mass of each protostellar core, depend- 
ing on its initial mass and the proper ties of the surrounding 
turbulent flow (|Bonnell et alJIioOlbl ). On the other hand, 
winds from the most massive stars play a crucial role for the 
energetic balance in the ISM of the cluster. Numerical sim- 
ulations have also confirmed that radiation feedback from 
young protostars can also help shape the cluster IMF, by 
setting a lower limits to the Jeans mass in the ambient gas 
(e.g. Offner et al. 2009). Effects of accretion and energetic 
feedback fr om winds are imp lemented in a model of the clus- 
ter IMF of iDib et all (j2010l ). 

To sum up, current theories of the IMF describe 
the latter as the combined result of gravoturbulent 
MC fragmentation and competitive accretion, offering 

different scenarios, responsi ble for its substellar part 

l|Bonnell. Larson fc Z inncckcr 120071 ). In this Paper we 
present a semi-analytical model that implements these two 
basic physical mechanisms and is in agreement with the gen- 
eralised multip art power-law form of the observational IMF 
(|Kroupa|[200ll ). The 'backbone' of our model is the mass 
distribution of gravitationally unstable, star-forming clumps 
which is linked to the probability density function (pdf) of 
the density field at different scales in MCs. 



1.1 Influential (semi-)analytical works on the IMF 

The most successfu l (sem i- ) analytical IMF models of 
Padoan fc Nordlundl ([20021') (hereafter: PN02) and of 
Hennebelle fc Chabrierl (|2008l ) (hereafter: HC08) are based 
on the pdf of the density field, resulting from supersonic tur- 
bulence. The distribution of the density n per unit volume 
is lognormal and the position of the maximum is a function 
of the standard deviation (stddev) a: 
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similarity at different scales ( N C (L) oc L ). More recently, 
IDib. Kim fc Shadmehril (|2007t ) extended the model of PN02, 
including dynamical evolution of the prestellar cores and 
their coalescence and collapse. The main criticisms to PN02 
model concern the imp lementation of th e geometry of the 
post-shock structures (|Elmegreenl l2007h and of turbulent 
and magnetic support in the collapse conditions (HC08). 
The mass of a formed dense core should also influence the 
suggested Jeans mass distribution. 

The extensive analytical model of HC08 is an attempt to 
refine the ph ysical description of collap sing structures in the 
framework of iPress fc Schechterl (1 19741 ) . A density threshold 
is imposed for collapsing clumps while turbulent and mag- 
netic support are implemented in the threshold conditions 
- through dependence of the turbulent Jeans mass from the 
mean turbulent velocity (V 2 ^) 1 ^ 2 an d a modification of the 
Mach number in the MHD case, respectively. Scale depen- 
dence of the pdf width is introduced in the natural lim- 
its between injection and dissipation in a turbulent cascade 
process. Accretion onto the dense prestellar cores is assumed 
to happen only within the initial gravitationally unstable re- 
gions and not from external sources. That arises consistently 
from the picture of fragmentation in collapsing structures 
with power-law density profile p oc r~ a (a < 2). No further 
fragmentation during the core collapse is assumed. 

Both mentioned models reproduce the high-mass, 
power-law part of the IMF with Salpeter slope F ~ —1.35 
from the implied turbulent theory. In PN02 this prediction is 
recovered from scale-dependence of K-ms and self-similarity 
at different scales, while in HC08 it is a substantial effect of 
inclusion of turbulent support. The IMF peak mass M v in 
PN02 model is a function of the average Bonnor-Ebert mass 
and scales inversely with the rms Mach number (or, the 
rms Alfvenic Mach number in the MHD treatment). Adopt- 
ing typical values of the ISM parameters in star-forming 
core^J, one gets M v ~ 0.16 — 0.25 M©. In the approach 
of HC08, the peak of the IMF is a result of transition be- 
tween the regime with significant turbulent support and the 
purely thermal behaviour. Its value is approximately con- 
stant (M p ~ 0.1 Mq) in models with insignificant turbu- 
lent support (small local Mach numbers) and varies within 
a larger range (~ 0.03 — 0.3 Mq) and toward lower masses 
in models with increasing supersonic turbulence. 
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where z = In (n/no) and no is mean density in the considered 
volume. The std dev depends on the Mach numbe r A4 for 
supersonic flows (|Padoan. Nordlund fc Joneslll997T ): 



a 2 = ln(l + b 2 M 2 ), 



0.5 



(3) 



In the semi-analytical model of PN02, the pdf of the 
density field is transformed into a distribution of Jeans 
masses. This approach serves for the probabilistic evalua- 
tion of the local collapse conditions for dense cores and hence 
to set up a low-mass cut-off of the IMF (see formula (24) 
in PN02). The cores themselves are identified as fragments 
of layers or filaments, formed by shocks in the supersonic 
turbulent flow. Their mass distribution is derived from the 
MHD shock jump conditions, whereas their number is con- 
trolled by a natural assumption for turbulent flows: self- 



1.2 Main features of the presented model 

Our semi-analytical model attempts to combine some of the 
mentioned ideas about the origin of the IMF with several 
refined physical assumptions about SF process: 

• Starting point is a turbulent cascade process in the ISM, 
which for simplicity we assume to be isothermal at T = 10 K 
at all scales. It leads to the formation of clumps at different 
scales L, beginning at an injection scale Li n j = 100 pc and 
proceeding downwards to Lq > 0.1 pc. Mean density and 
velocity dispersion at each scale are estimated according to 
the Larson's (1981) relationships. 



1 Mean column density 
10 K, magnetic fiel d B 
iPadoan etail ||2007|) . 



10 22 cm" 2 , n ~ 10 4 cm~ 3 , T = 
10 fiG; see formulae 9) and 10) in 
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• The clumps, generated through the supersonic turbu- 
lent flow at a given scale, have a lognormal density distribu- 
tion that corresponds to the pdf of the density field. Assum- 
ing a power-law relationship between their mass and size, 
the clump mass distribution is obtained. 

• Gravitationally unstable clumps are selected at scales 
La ^ L ^ Li n j and their composite clump mass function 
(CMF) is derived. Each of them fragments further and pro- 
duces prestellar cores at the local Jeans scale. 

• Competitive accretion onto the formed prestellar cores 
distributes the mass of the remaining gas from the frag- 
mented clump. Clumps, which have transformed all their 
material to prestellar cores (by fragmentation and accre- 
tion), are replenished at constant rate. 

• Eventually, the IMF is derived from the compos- 
ite CMF, assuming ongoing generation (replenishment) of 
clumps and taking into account their fragmentation and the 
accretion on the formed prestellar cores. 

In Section 2 we describe how the composite mass func- 
tion of gravitationally unstable clumps is derived. Section 3 
presents the derivation of the IMF from the CMF through 
fragmentation and accretion. We discuss the model predic- 
tions in Section 4 and summarise the main contributions of 
our model in Section 5. 



2 THE CLUMP MASS FUNCTION (CMF) 

2.1 The concept of clump and its physical 
parameters 

Supersonic isothermal turbulence in MCs creates an in- 
tricate network of interacting shocks resulting in den- 
sity fluctuations, described statistically by a lognormal 
pdf (Eq. [TJ. This is to be expected from the stochastical 
nature of turbulent flows and is confirmed from various 



numerical simulations (IPadoan. Nordlund & Jonesl 19971; 


Ostriker, Stone & Gammie 


200ll; Li. Klessen & Mac Low 


2003 


; Kritsuk et al. 12007 


; IFederrath. Klessen & Schmidt 


2008 


; iFederrath et al.l 20101) . To attribute a given densitv 



from the pdf to a spatial physical object (condensation) 
with particular volume and shape, one needs an identi- 
fication scheme for stru ctures in MCs. Common meth- 
ods like CLUMPFIND (|Williams. de Geus fe Blit3 Il994f ) 
identify contiguous structures in datacubes with densi- 
ties over a given threshold as distinct condensations, la- 
beled generally "clumps". Usually, clumps are a broadly 
defined group of objects with sizes from several tenths 



of parse* 
10 



t o 5, 10 P c and with masses in the range 



10 4 M m (iKramer etal.ll 19981 ; I Williams. Blitz fc McKeel 



|2000| ; iKauffmann et alj l2010af T They are, at least, con 
fined by the external p ressure of the ambient medium 
(|Bertoldi fc McKedfl99a ) and those which are virialised or 
gravitationally bound are the massive clumps where cluster 
formation takes place (|Williams. Blitz fc Starkll 19951) . 

Hereafter in this Paper, we use the term clumps for 
MC fragments of various shapes that have been formed by 
turbulent shocks at scales 0.1 < L < 100 pc and have scale- 
dependent density distributions according to Eq.Q] They are 

2 In order to be distinguished from dense prestellar cores with 
sizes < 0.1 pc. 



potential SF sites where clusters of dense prestellar cores 
could form (Sect. 13 - 1 P . A statistical approach is adopted, 
assuming for simplicity that all clumps have cubic shape - a 
clump size I corresponds to clump volume of I 3 . We suppose 
a power-law mass-size (m — I) relationship for clumps which 
also implies a mass-density (n — m) relationship: 



/ m 



Ox 
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m 
m 



I \ 3/(1-*) 



(4) 



where no(L), mo and lo are units of normalization. The 
existence of a mass-density (or, mass-size) relationship for 
molecular clouds and clumps can be expressed in terms 
of the combination of a velocity scaling law with differ- 
ent c onsiderations of energy balance: equipartition o f ener- 
gies (|Ballesteros-Paredes fc Vazquez-Semadenilll99 5), virial 
equilibrium or relationships between the virial parameter 
or the Jeans numb er and the clump mass (|Dib et al.ll2007l ; 
IShettv et al.ll201(il ). Depending on the chosen approach, the 
result for the power exponent is in the range —1.5 < x < 0.4 
(its demonstration is beyond the scope of this paper). Rela- 
tively large variations of x are obtained as well from present 
observational studies of cloud fragments in MCs (—2 < x < 
—0.5, Kauffmann et al. 2010b) and numerical simulations 
(-2 < x < 0., Shetty et al. 2010). Therefore x is taken to be 
a free parameter of our model. To avoid implementation of 
additional physics, scale independence of x is also assumed 
although this might be a crude approximation. 



2.2 CMF at given scale 

2.2.1 Clump density distribution 

As mentioned above (Sect.[T]2}, we use the proposed clump 
mass-density relation (Eq. [4| to derive the CMF from the 
clump density distribution. An appropriate statistical de- 
scription of the latter is a lognormal pdf (Eq. [1} . Its scale 
dependence is set by the scaling of the Mach number (cf. 
Eq. [3J, i.e. of the velocity dispersion Kms(L), and by the 
choice of normalization unit no(L). Observational estimates 
of those quantities for the large range of considere d scales 
could be provided from the so called "Larson's laws" (|Larsonl 
[l98l]). 



V„ 
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1 pc 



[km / s 



0.4, 



n = (n) = 0.34 [10 4 cm 

1 PC/ 



(6) 



(7) 



These relations do not imply anything about the specific 
parameters of the clumps themselves and their behaviour 
(e.g. the mass-density relationship) but are used here to give 
appropriate, statistically averaged values of V Trns and no at 
each scale of the turbulent cascade. Choosing no = (n) corre- 
sponds to setting a typical density in the considered volume, 
so that clumps could be thought as contiguous structures 
within isodensity contours in units no. 

The other factor that determines the clump density 
distribution is the parameter b (cf. Eq. [5). Its value is 
subject to some uncertainty. Observations (|Bruntl I2010T) 
and 3D numerical simulations (|Padoan. Nordlund fc Jonesl 
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[1993; 



Kritsuk et all 120071) suggest . 25 ^ b ~ 0.6. How 



Federrath. Klessen fc Schmidtl (|2008T l demonstrated 



that the type of turbulence forcing should be taken into 
account when b is evaluated - compressive forcing produces 
a three times larger stddev of the pdf than in the purely 
solenoidal case (see also Federrath et al. 2009). That corre- 
sponds to variations of b in the range 0.33 — 1.0 which we 
adopt further in this work. 



2.2.2 Clump mass and size distributions 

After the scaling and the normalization of the lognormal 
clump density distribution have been specified, one is able 
to derive the mass and size distributions. Those are also log- 
normal, as Eq.|4]imply, and we denote them p m (ln(m/mo)) 
and pi (In (l/lo)), respectively. Their parameters are obtained 
straightforwardly: 

z m = \n(m/m ) : z m , max = Z max /x , o m =cr/\x\ (8) 

Zl=\n{l/h) \ 2!, max = Zmax/Ca; , CT( = <t/\C x \ (9) 

A lognormal pdf of the density field, generated by tur- 
bulence, has no natural limits; thus the derived clump mass 
and size distributions have not such. There is no natural 
choice for the normalization units mo and Iq as well. We ob- 
tain the size and mass range and the normalization units by 
use of the requirements for volume and mass conservation: 



L 3 = N^2l 3 (Ni/N) ~ Nil [ e 3z; 



p{zi)dzi (10) 



{p)L 3 = [in L 3 = M sca i ~ Nm 3 I e Zm p(z m ) dz„ 



(11) 



where Ni is the number of clumps of size I, N is the total 
number of clumps in the volume V and the mean molecular 
weight is taken to be u = 2Am p . Consecutive numerical 
integration, performed symmetrically around z;, m ax, yields 
the size limits of clumps as lower and upper cutoffs when 
Eq. (|10[) and (|11|) are satisfied. The integrals on the right- 
hand sides have analytical solutions in the limits (— oo, +00). 
Using them, one obtains an approximate relation between 
the normalization units: 



m 



exp 



1 - x 



(12) 



For the sake of numerical integration, it is appropriate to 
choose Iq to be small and proportional to scale size L (of 
order of few percent of it): 

lo = hL 



2.2.3 Mass distribution of unstable clumps 

Clumps for which the mass exceeds the local Jeans mass 
mj(n) are, by definition, gravitationally unstable. Their 
fraction or formation at given scale L depends on some crit- 
ical mass m CI (L) that separates the regimes m/mj(n) = 
m/mj(m) ^ 1 and m/mj(m) < 1. 

Introducing the chosen normalization unit no = (n) into 
the Jeans mass expression mj ~ 1.5 (T/W K) (n/10 4 cm" 3 ) 
and using the mass-density relationship (Eq. [4]), we obtain 
the criterion for gravitational instability of a clump: 



K(x,L) 
mj yl pc ; 

which yields a critical mass: 



cr = K{x,L) *+2 

1 1 pc 



m 



> 1 



(13) 



(14) 



The factor K(x, L) varies within an order of magnitude for 
x > and about 2 orders of magnitude for x ^ — 3. A brief 
analysis of Eq.[l3]and[T4]displays two special physical cases. 

• x = —0.5: That is the Larson-like case in which n oc Z _1 
(cf. Eq.HandEJ). The critical mass is huge (~ 700-1500 M ) 
and depends weakly on the scale: 



= K{x,L)s 



(15) 



* x = —2: All formed clumps are unstable by definition 
(m oc 7i~ 1//2 ) and an appropriate name is Jeans-like case. 
It is obviously unrealistic, although the mass range of un- 
stable clumps is significantly larger than in the Larson-like 
case. The criterion for gravitationally instability (Eq. 1 1 3 p 
is fulfilled only for scales where K(x,L) > (L/l pc) _3//2 , 
or L > 70 pc. If realised (partially) in nature, the Jeans-like 
case would correspond to starburst regions at larger galactic 
scales. 

The clump mass distributions at different scales, pro- 
duced in the above cases, are shown in Figure[T](left panels). 
Two other cases worth mentioning are: 

* x — 0: Homogeneous medium with no density fluctua- 
tions (case of decayed turbulence). 

* x ~ 1: Degeneration in the clump size - clump density 
relation (Eq. [3]), i.e. clumps of approximately constant size 
have densities in an extremely large range. 

The clump mass distributions for some other values of x 
are plotted in Fig. [1] (centre and right panels) . Apparently, 
turbulent cascade in clouds of sizes 1 — 100 pc gives birth to 
gravitationally unstable clumps in a large range of masses 
1 — 10 5 Mq. A lower mass limit of < 1 Mq is achieved for 
positive x < 0.5. Negative values of x below the Jeans-like 
value (< —2.) condition that almost all clumps that would 
form at large scales [L > 80 pc) are unstable. That leads to 
a top-heavy composite CMF we derive in the next Section. 

The abundance of gravitationally unstable clumps de- 
pends not only on the scale L (Eq. I14p . but also on the 
mass distribution parameters: the position of the maximum 
and the width (Eq. [2] [3] and [Sj . The latter dependence be- 
comes stronger as \x\ grows - Zmax,m at a given scale shifts 
to lower masses while a m decreases (see Fig. [TJ|. In such 
cases and especially at small scales, shape and position of 
the mass distribution are the crucial factor that determines 
whether the formation of unstable clumps would occur or 
not. In other words, the formation and the fraction of un- 
stable clumps at given scale L are a product of interplay 
between the values of x, m cr and a m . This is illustrated in 
Fig. [2] Unstable clumps form in the mass range: 

m CI (L) ^ m > rn°™ tof! (x,L), x < -2, (16) 

m CI (L) s£ m< m^ toS (x,L), x > -2, (17) 

where mass cutoffs (mJ.°™ off , rn^ toB ) are determined from 
simultaneous numerical integration of Eq. [10] and [11] 
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Figure 1. Clump mass distributions for different values of x and at different scales: 100 pc (red), 70 pc (green), 50 pc (blue), 30 pc 
(grey), 10 pc (violet) and 1 pc (black). The forcing parameter is fixed: b = 0.5. The portions that correspond to gravitationally unstable 
clumps are shown with bullets. 



(Sect. 1223) . If m cr (L) < m}°^ oS (x,L) (or m cr (L) > 
Tn^ toS (x, L)) for given pair (x,L), no unstable clumps are 
formed. As seen in Fig. [2 small positive values of x yield 
lower scale limits of formation of unstable clumps Lo(x): 
Lq ~ 1 pc for x ~ 0.5 and a value of Lq ~ 0.1 pc, typical 
for transition from clumps to dense cores, is achieved only 
for < x < 0.25. On the other hand, negative values of x 
lead to formation of unstable clumps only at scales > 10 pc. 
Noteably, turbulent flows remain supersonic at all scales, for 
all choices of x. 



2.3 Composite CMF 

The obtained mass distributions of unstable clumps have to 
be summed up over the range of scales Lq(x) < L < L- m j to 
derive a composite CMF, representative for star formation 
at galactic scales. In fact, aggregates of MCs of different sizes 
are embedded in the general ISM. If all of them are gener- 
ated by supersonic turbulence through a cascade process, 
so a self-similar distribution in a given volume unit can be 
assumed. Hence the total number of clumps TV oc L~ 3 (see 
PN02) and the contribution of clumps, generated at scale L, 
to a logarithmic mass bin d(lgm), scales as 



N(L, d(lgm)) = ( 



jV(L inj , d(lgm)) 



(18) 



where Li n j is the chosen injection scale of the turbulent cas- 
cade. Then, the total number of clumps, contributing to a 
mass bin, is: 

JN/tot d(lg m) = J N(L, d(lgm))dL 

L (x) 



The derived composite CMFs for a set of exponents x axe 
plotted in Fig. [3] Two different types of behaviour are appar- 
ent. For negative values of x, implying a steeper density-size 
relation than in "Larson's second law" (n k 1°, a < — 1), 
the CMFs are top-heavy, asymmetric, with steep edges 
and shallow power-law parts. The slopes are in the range 
—0.6 ^ r ^ —0.8 , in agreement with observations of MC 
clumps, associated with active SF regions (e.g. Kramer et 
al. 1998), but the turn-over mass M to , cmf is significantly 
larger, by an order of magnitude. The cases x > —2 yield 
smooth, power-law CMFs spanning a wide range of masses 
and with slope T « —1.2, close to the Salpeter value. The 
turn-over mass is < 1 Mq for small positive values of x 
and depends on the low-mass limit of the CMFs. The latter 
is determined by the lowest scale of unstable clump forma- 
tion Lo(x). Therefore, in our approach, Mto, cmf is a scale- 
averaged, model-dependent (on the value of x) and not an 
universal phenomenon. 



3 FROM THE CMF TO THE IMF 

The derivation of the IMF from the CMF reflects the pro- 
cesses of fragmentation and further accretion on the gen- 
erated protostellar cores. Hereafter, we will use the term 
prestellar cores (or, simply, cores) to denote those objects 
that form via the fragmentation of unstable clumps, and 
denote their density and mass with n' and m', respectively. 



3.1 Fragmentation of the unstable clumps 

The most straightforward approach to fragmentation is to 
assume that each unstable clump of mass m fragments into 
(m/mj) protostellar cores. That procedure, however, is not 



6 Veltchev, Klessen & Clark 



10 H 



= 1 10 2 



a o 

E 



10 u 



10 



-2 



Unstable clumps 













1(T 



10" 



10 



o 

U -j o 



10 



-2 



-6 -4-2 2 

Exponent x 

Figure 2. Mass ranges of gravitationally unstable clumps in the 
regimes separated by the Jeans-like case: x < —2 (left) and x > 
—2 (right). Unstable clumps form if m cr > "i'°^ tof j for x < —2 
and «i"u ttof j > "icr for ie > —2.. The scales are designated like in 

Fig. El 



well justified physically due to, at least, two reasons: (1) 
it implies unrealistically high SF efficiency; (2) it neglects 
the internal structure and dynamics of a collapsing clump. 
I Goodwin et al.l (|2008l ) adopt a purely probabilistic approach 
to describe fragmentation of low mass clumps (0.1 — 100 Mq) 
into 2 or 3 fragments, assuming constant ratios of binaries 
to triples and binaries to single stars, constant SF efficiency 
and an uniform probability distribution. We consider clumps 
in a significantly larger mass range and, in most cases, with 
(m/mj) 3> 1. Therefore we take a different approach on 
fragmentation. 

In view of the fractal nature of turbulence, it is reason- 
able to describe clump fragmentation (also) through a log- 
normal core mass distribution. This should be a distribution 
of Jeans masses, i.e. local density contrasts n'/n(raj) within 
the clump generate colapsing cores with corresponding Jeans 
masses m' = mj oc (n') _1//2 . The peak of the distribution is 
assumed to be at the clump Jeans mass mj(n): 



mj + Am(ff) 



j m'j N(m'j) dm'j ; N(m' 3 ) = 



771 j — A m { cr ) 



p(mj + Am(a)) 
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Figure 3. Composite CMFs, derived for x ^ —2. (top, open 
symbols) and x > —2. (bottom, filled symbols) values of x, like 
in Fig. [T] —2.0 (Jeans-like case; circles), —3.0 (triangles), —6.0 
(diamonds); —0.5 (Larson-like case; pentagons), 0.25 (squares), 
0.5 (squares), 2.0 (circles). The CMFs in the bottom panel are 
artificially shifted vertically for cla rity. Typical range of slopes 
(grey area) for observational CMFs jKramer et alj|l998h and the 
best-fit slope for x ^ —0.5 (dashed line) are plotted. 



where Am(a) gives the numerical integration cutoffs. The 
width of the cores' mass distribution a m i = a n i /2 (x = — 2) 
is calculated from Eq. [3] In that equation, the global-scale 
Mach number M(L) has to be replaced by the local Mach 
number Mi, reflecting the velocities w rms at scales of the 
clump size I (see comments on that in iBonnell fc Bate! 
(|200q) . Sect. 2). Clump sizes in our model span the range 
0-01 < I < 10 pc which includes a variety of different phys- 
ical regimes. Turbulence dominates the observed linewidth 
down to sizes of I ~ 0.1 — 0.2 pc |Barranco fc Goodman! 
Il998h . Its scaling in that regim e obeys relation: f rm3 ~ l a , 
with a slope a ~ 0.41 - 0.43 (jPadoan et al.1 l200fj , 120091 ') . 
close to Larson's value (Eq.O, or steeper: 0.48 ^ a ^ 0.75 
(|Hever. Williams fc BrurraboOfJ l. At sizes 0.1 - 0.2 pc the 
thermal component of w r ms becomes comparable to the tur- 
bulent one and the slope a decreases to 0.1 — 0.15. Even- 
tually, at scales ~ 0.04 pc, "coherent co res" with w rms ~ 
const are observed (|Goodman et al.lll998h . An appropriate 
parametrization of the velocity dispersion in the fragmenting 
clumps, describing a smooth transition between turbulent 
and thermal regimes, is: 



M = vZ(c + (i/i ) 2a ) , 



(19) 



where we take vq — 0.8 and a = 0.48, obtained fo r MC 
with low SF efficiency (|Hever. Williams fc BruntJl2006l ) . and 
fiducial size lo = 0.8 pc. C is a small constant for which we 

choose a value of 0.1. 

Similar to I Goodwin" et al.1 (|2008l ), we define prestellar 
core formation efficiency e (PCFE) as the percentage of 
clump mass m that remains bound in the cores, after the 
initial fragmentation. According to observational estimates, 
the total gas mass in yo ung stellar c lusters is typically many 
times the mass in stars (|Ladalll99ll ). Therefore it is reason- 
able to choose small values of e (see also Krumholz & Tan 
2007, ApJ). To simplify the calculations, we adopt constant 
values of this par ameter, although one can expec t that it is 
scale-dependent (|Vazquez-Semadeni et all [20031 ) . The rest 
of the initial clump mass m gas (to) = (1 — e)m is distributed 
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further among the formed cores through competitive accre- 
tion. 



3.2 Description of accretion on the prestellar 
cores 

Accretion in protostellar clusters proceeds in two main 
phases: gas-dominated poten tials and stellar-dominated po- 
tentials (jBonnell et al.|[2001bl ). In our approach, small values 
of e imply that the gas dominates the potential within the 
fragmented clump at some initial moment to- Description of 
accretion on cores under such conditions is difficult to imple- 
ment in a (semi-) analytical approach. The general formula 
for the accretion rate on a protostar of mass M is: 



M ~ irpv Te iR a 



(20) 



where p is the gas density, v Ic i is the relative gas-protostar 
velocity and i? a cc is the accretion radius. If gas dominates 
the potential and the collapse is isothermal , a good esti- 
mation of 7? acc is the tidal-lobe radius JJtidai (jBonnell et alj 
12001 j ). However, the evaluation of .Rtidal = f(r) for each 
core requires knowledge of its current position r within the 
clump, i.e. further assumptions about cores' spatial distri- 
bution and its evolution in time. Also, tidal accretion tends 
to dominate in the fragmentation phase (as described in 
the previous Subsection ), not in the N-body dynamic s that 
takes place afterwards (|Bonnell. Clark fc Bate! [20081 ) . It is 
therefore more a ppropriate to use the Bondi-Hoyle accre- 
tion description (jBondi fc Hovldfl944h , devised for stellar- 
dominated potential: 



Ra 



Rbh 



2GM 



M « Me 



47rp 



(GM) 2 



K 2 ol + cf)3/2 



(21) 



(22) 



where the gas density in protostar's vicinity is assumed to 
be uniform. The latter assumption is not far from reality 
also in a moving gas medium, i f the gas distributi on re- 
tains its general form with time (jBonnell et al.l[l997l ). The 
density profile p(r,t) within a collapsing clump approaches 
quickly the isothermal for m (p oc r~ 2 ), with a small nucleus 
of nearly uniform density (jLarso n 1969). One can adopt an 
approximation of the spherically averaged gas density profile 
in fragmented and collapsing clumps: 



P(r,t) 



Pc(t) 

Pc(t) (r B (t)/r 



r < r c (t) 
r > r c {t) 



(23) 



where r c (t) is the radius of the clump nucleus with uniform 
density p c (t) at a fixed moment of time. The time evolutio n 
of the latter quantity is approximately linear (jLarsonlll969h : 



p c (t) = p c (t ) + At 



(24) 



whereas r c (t) could be obtained from the equation of gas 
mass m gas (t) = f Aivp(r,t)r 2 dr and hence depends on the 
accretion rates within the clump. Then the mean density 
within a clump of size / and at a fixed moment t is: 



Pit) = \ J P (r,t)dr = 2 Po (t) r -^(l 



Tc(t) 



(25) 



This estimate is physically more correct than to assume 
decreasing (due to accretion) uniform mean density. The 
relative gas-to-protostar velocity could be derived from the 
clump Mach number Mi (see previous Subsection). Eventu- 
ally, accretion on a protostellar core of mass m! is described 
in our approach through the Bondi-Hoyle formula: 



m'(t) = 4irp(t) 



(Gm' 



(« r 2 el (M) + ci) 3 / 2 



(26) 



and the total mass of the gas within the clump is being 



exhausted at a rate m E 



fm'N(m')dm'. 



3.3 Derivation of the composite IMF 

3.3.1 Parameters and timescales 

The derived IMF depends on three basic parameters. Two 
of them determine the CMF: the exponent x in the clump 
mass-density relationship and the turbulent forcing param- 
eter b (Eq. |3J; while the third, the PCFE e, acounts for 
fragmentation. 

The free-fall time rg of a unstable clump is a natural 
measure for description of its evolution (fragmentation and 
internal accretion). Since 



Tg 



\32Gp J 



1/2 



-1/2 -x/2 

oc n ' oc m ' 



the range of clump evolution timescales is determined by 
the clump mass distribution and, hence, from parameters x 
(Eq. [8j and b (Eq.[3]). The maximal value Tg t max is achieved 
at the largest (injection) scale L- 1Yl - s — 100 pc and it could be 
taken as a characteristic evolution time of protostellar cores 
generated through the turbulent cascade. The dependencies 
of Tg on clump mass m, turbulent scale L and on x are 
illustrated in Fig. [4] Larger scales produce wider spans of 
free-fall times and this correlation is very sensitive to the 
chosen value of x. For |x| > 1, unstable clumps that are 
produced at scales L > 50 pc would evolve at free-fall times 
varying within two orders of magnitude! Moreover, Tg (m, L) 
are comparable to or even greater (in the Larson-like case) 
than the turbulent crossing time at the given scale t ci (L) oc 
L 0,6 . That makes models with such values of x implausible 
for description of the SF process. On the other hand, the 
cases 0.2 < x < 0.5 produce a range of clump free-fall times 
from few tenths to several Myr, which is consistent with the 
rapid SF model, in its observational or theoretical timescales 
(Ballesteros-Paredes, Hartmann & Vazquez-Semadeni 1999; 
Pringle, Allen & Lubow 2001). 

The observational IMF at galactic scales is represen- 
tative of stars born at different scales and hence within a 
large range of clump free-fall times Tg (n, L). In some of the 
protostellar clusters the gas is exhausted and hence the ac- 
cretion is halted. Others are still embedded in their natal 
cloud which sets certain observational limits to their de- 
tection and/or estimation of their masses. In this study, we 
abstain from implementing such observational aspects in the 
model and take into account only fragmented clumps (pro- 
toclusters) in our derivation of the IMF with initial mass 
em where the accretion has halted due to gas exhaustion at 
a moment of time t aC c, less or comparable to Tg. The SF 
efficiency (SFE) at fixed time t is defined as the ratio of the 
total (initial + accreted) mass of the prestellar cores M* (t) 
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Figure 4. Free-fall times of unstable clumps produced at scales 
100, 70, 50, 30, 10 and 1 pc (designated as in Figs. 1 and 2) 
and for different values of x. Note the weak scale dependence for 
x'| < 0.5. The range of turbulent crossing times r cr (L) is drawn 
for comparison (yellow area). 



Table 1. Maximal clump free-fall times and SFEs. 



X 


b 


'ff, max 

[Myr] 


SFE at t = T a , 
e = 0.10 t = 0.40 


max 

t = 0.70 


0.25 


0.33 


5.4 


0.12 


0.51 


0.89 


0.50 


0.33 


7.5 


0.56 


0.75 


0.89 


0.97 


0.33 


13.4 


0.20 


0.36 


0.60 


0.25 


1.00 


5.0 


0.12 


0.66 


0.90 


0.50 


1.00 


5.8 


0.30 


0.55 


0.79 


0.97 


1.00 


15.0 


0.20 


0.35 


0.58 



to the total mass M to t involved in the whole turbulent cas- 
cade. The SFEs, obtained for a set of values of x and e and 
at t — Tff imax , are given in Table Q] 



3.3.2 Treatment of the timescale problem 

To account for the various accretion timescales i acc within 
the clumps i n derivation of the IMF, we follo w an approach 
discussed bv lClark. Klessen fc Bonnelll (|2007l ). Without any 
consideration of fragmentation and accretion, those authors 
demonstrate that the IMF would differ substantially in form 
from the progenitor CMF (existing over timescale tcmf), if 
the evolution timescale for the clump i ovo i depends on its 
mass m: /imf ~ /cmf X tcmf /t cvo i(m). In our model, stel- 
lar progenitors are the cores within a fragmented clump, 
with protostellar core mass function (PCMF) /pcmf(ctti) 
and evolution timescale i acc (em). Permitting the clump pop- 
ulation to be constantly replenished, we have for the local 
IMF (LIMF) at t > t acc (em), produced from the PCMF 
through competitive accretion: 



/limf , L) 



tacc(em) 



/pcMF(em,L) 



(27) 



Eventually, the composite IMF at a fixed moment of 
time t is derived through integration over the range of scales 
where unstable clumps are formed: 



/ 



/limf , L) dL 



(28) 



L (x) 



4 RESULTS AND DISCUSSION 

Plausible evolution timescales of the protostellar clumps 
(~ Tff(x)) are obtained for the slopes of the clump mass- 
density relation (Eq.[4j) in the range < x < 1 (cf. Table [TJ. 
Therefore we focus on the results for the IMF in the cases 
x = 0.25, x = 0.5 and x = 0.97 « 1), plotted in Figs. [51171 
Several apparent features should be pointed out: 

• A Salpeter slope of the high-mass IMF is generally re- 
produced for small PCFE e for all choices of the exponent 
x. This result is not trivial since it is not a direct effect 
of transforming the CMF with a similar slope (cf. Fig. |3J), 
but includes the effects of fragmentation and accretion. On 
the other hand, steeper slopes are obtained as e is increasing. 
That is understandable, because decreasing of the initial gas 
mass m gas in the clumps minimises the role of accretion on 
the protostellar cores in reshaping their mass distribution. 

• Intermediate-mass IMF with shallow (negative) slope 
and a mass range within an order of magnitude (0.06 < 
M < 0.6 Mq) is derived mainly for solenoidal turbulent 
forcing. Compressive forcing (b = 1.0) tends to produce a 
narrow IMF, with a peak shifted toward the BD mass range. 
The apparent exception is the case x — 0.25: both extreme 
regimes of forcing yield an IMF in a good agreement with 
observations. 

• A discontinuity of the IMF is evident in the BD mass 
range for x = 0.50, 0.97 and higher PCFE (Fig. [6] and [7J. 
It is about o ne order of magn i tude o r less higher than what 
is found by iThies fc Kroupal l|2007h in the young clusters 
Trapezium and IC 348. 

The model with x = 0.25 and low PCFE exhibits the 
best consistency with the observational IMF. According 
to Eq. 0, this case corresponds to an interesting clump 
geometry - its density is exactly proportional to the lin- 
ear size: n ex. I. This result is in apparent discrepancy 
with "Larson's second law" (Eq. [7}. It should be consid- 
ered cautiously but not automatically rejected. The rela- 
tion n oc i -1 (or, equivalently, m oc l 2 ) is representative 
for clumps with column density constant within an order 
of magnitude and thus seems to be ob servationally biased 
(IBallesteros-Paredes fc Mac Lowll2002l ). Also, it holds for 
clumps that are in equipartition between the gravitational 
and the kinetic energ y or between the gravita tional and 
the magnetic energy (|Ballesteros-Paredes! |2006l ) . Different 
considerations of the clump energy budget yield a variety 
of mass-density relationships, with —1.5 < x < 0.4, as 
we intend to demonstrate in a forthcoming paper. More- 
over, from the perspective of observational clump identifi- 
cation schemes, the Larson-like case m oc I 2 (x = -0.5) 
sets a n upper limit on slopes of clump mass-size relation- 
ships jKauffmann et alj|2010a) while analysis of 3D clumps 
from numerical simulations ( Shettv et alj [20 10;) increases 
this limit up to x > 0, depending on the size range. 

The obtained characteristic evolution timescales and 
SFEs (Table [TJ for the model x = 0.25 (low PCFE) seem 
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Figure 5. Initial mass function for x = 0.25. Results for purely 
solcnoidal forcing (left panels) and for purely compressive forcing 
(right panels) are plotted: a) (top) fixed, low PCFE e = 0.1 and 
time evolution 0.2, 1.0, 1.5 Tg max ; b) varying PCFE (bottom): 
e = 0.1 (red), e = 0.4 (blue), e = 0.7 (violet) , at t = Tg max- 
Multipart power-law IMF according to iKroupal j200ll) is plotted 
for comparison. 
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Figure 6. The same like Fig. \5\ but for x = 0.50. The range o f 
slopes of the BD-like IMF as obtained bv lThies fc Kroupal j2007l l 
is drawn as grey area. 



also plausible. In that case, free-fall times of the unsta- 
ble clumps are about an order of magnitude less than the 
turbulent crossing times r cr (L) (Fig. The SFEs are in 
agree ment with estima tes of about 5-10 percent in giant 
MCs (|Clark et al.ll2005h , although the latter could vary in a 
wider range, depending on the d ynamical state of the cloud 
l|Clark. Bonnell fc Kless"enll200Sl ). Accretion rate calculated 



Figure 7. The same like Fig. [6] but for x = 0.97. 

according to the chosen description (Sect. 13. 2[) also fits the 
current state of the theory. Nu merical simula t ions of accre- 
tion in protostellar clusters by iBonnelT et all (|2001al ) show 
that most of the gas is exhausted at t acc ~ Tg. This sets 
up only a lower limit of t acc since feedback from massive 
stars and magnetic fields were neglected in their work but 
is an appropriate reference value in view of the description 
of accretion we adopted. In our model, the majority of the 
unstable clumps at all scales that have given birth to larger 
clusters (> 20 cores) have consumed the remainining gas at 
t — Tff )Inax . Thus the shape of the IMF is practically fin- 
ished at timescales about t = Tff. max (Fig. [5] top pannels) 
and further accretion does not influence it significantly. 

An apparently realistic IMF is derived from models 
with x = 0.50 and x = 0.97 and higher PCFE (Fig. [6] 
and bottom panels). However, such values of the expo- 
nent x cannot be substantiated theoretically and lead to 
clump mass-size relationships in drastic contradiction both 
with present observational studies and numerical simula- 
tions (|Kauffmann et al.l [201Qbl ; IShettv et al.l |2010h . More- 
over, these models fail to reproduce the substellar IMF - 
discontinuity in that mass range is about an order of magni- 
tude larger than found bv lThies fc Kroupal (|2007h in young 
Galactic cluster^]. 

Accretion on protostellar cores plays a significant role 
for shaping the final IMF in our model. Therefore it is in- 
st ructive to make a brief comparison with the recent work 
of iDib et all (|2010l ) in which accretion is also one of the key 
factors and the IMF is also derived from a sum of local core 
mass distributions. Without going into details, we stress the 
impor tant conceptual differences. The model of IDib et all 
(|2010l ) is aimed at reproduction of the IMF in young clusters 
where the SF process has been rapid. The 'building blocks' 
in their approach are the prestellar core mass distributions 
derived through the formalism of PN02 for different loca- 
tions within a protocluster clump. Subsequent accretion is 



If such feature is going to be confirmed. 
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also described locally: on cores, injected in the protocluster 
at a uniform rate and in different epochs, until the gas in 
the protocluster is dispersed by stellar winds from massive 
stars. Our model of the IMF is aimed to be representative 
of star formation in galaxies, in a wider range of timescales 
(cf. Fig. |4}, and its 'building blocks' are unstable clumps 
generated through a turbulent cascade spanning a range of 
3 orders of magnitude in spatial scales. Their fragmentation 
into prestellar cores proceeds in a self-similar way, depend- 
ing on the local Jeans mass. The time evolution of the IMF 
in our model is based on the assumption that the clump 
population is constantl y repl enished (Sect. I3.3.2|l . Unlike 
the model of iDib et all (|2010l ). accretion is spatially aver- 
aged for the clump in consideration and feedback from newly 
formed massive stars is neglected. In view of the large clump 
statistics assumed in our model, such approach would not 
affect significantly the predicted IMF although refinements 
are necessary when the SF process is described specifically 
at the protocluster scale. 



5 SUMMARY 

In this work we developed a semi-analytical model of the 
IMF that takes into account the basic mechanisms in the 
SF process: gravoturbulent fragmentation and subsequent 
accretion on the prestellar cores. It is based on the idea 
of a turbulent cascade that starts from large injection scales 
~ 100 pc and transfers energy down to dissipation scales be- 
low 0.1 pc. Turbulence is essentially supersonic at all scales 
which results in network of interacting shocks and forma- 
tion of condensations (clumps) with lognormal distribution 
of densities. The main assumption in our model is the exis- 
tence of a power-law relation between the clump mass and 
the clump density n oc m x , where a; is a scale- free parame- 
ter. It leads to lognormal clump mass distribution and to a 
scale-dependent critical mass that determines the fraction of 
gravitationally unstable clumps (if any) at given scale. Sum- 
mation of their mass distributions over all scales yields their 
composite CMF: asymmetric with shallow power-law part 
(slope —0.8 < r < —0.6), for x < —0.5, and steeper, with 
a Salpeter-like slope (F ~ —1-2), for x ^ —0.5. The frag- 
mentation of unstable clumps is assumed to produce a log- 
normal mass distribution of prestellar cores, centred around 
the clump Jeans mass. Adopting constant core formation 
efficiency e, we let the rest of the clump mass (1 — t)m to 
be distributed on the formed cores through competitive ac- 
cretion. The duration of this process is within ~ 1 clump 
free-fall times rg(x) as the gas density in the protocluster 
approaches an r -2 profile. Eventually, the IMF is derived 
through summation over all clumps with halted accretion, 
build at all scales, and allowing constant replenishment of 
clumps with faster evolution. 

The free parameters of the proposed IMF model are the 
exponent x in the clump mass - clump density relation, the 
turbulent forcing parameter b and the core formation effi- 
ciency e. Models with x = 0.25 and low e turn out to be con- 
sistent with the observational IMF. The physics behind the 
proposed relation n ~ m x is to be substantiated. Fragmen- 
tation of unstable clumps is another element of th e model 
that n eeds further elaboration. As demonstrated bv lKlessenl 
(2001). it depends on the turbulent injection scale and hence 



may yield a variety of mass distributions of dense cores. 
Eventually, we point out that the implemented description 
of competitive accretion within a fragmented clump gives 
only general estimates. Accurate modelling of this complex 
process is practically impossible in an analytical or semi- 
analytical framework. 

To conclude, the proposed model reproduces correctly 
the effects of gravoturbulent fragmentation and competitive 
accretion on the IMF at galactic scales. It is based on reason- 
able physical assumptions about the SF process and should 
be considered as a first step toward a more extensive theo- 
retical framework. 
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